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I.  INTRODUCTION 

Our  previous  numerical  simulations  of  the  nonlinear  evolution  of 
the  collisional  Rayleigh-Taylor  instability  in  the  nighttime  equatorial 
ionosphere  (spread  F)  were  confined  to  small  (~  3km)  horizontal  scale 
length  initial  perturbations  and  hence  to  furly  developed  spread  F 
"bubbles"  of  approximately  the  same  size  in  horizontal  extent  [Scannapieco 
and  Ossakow,  1976;  Ossakow  et  al . ,  1979],  although  spatially  large 
vertically.  However,  observations  by  McClure  et  al  [1977]  also  indi¬ 
cate  ionospheric  ion  density  "biteouts"  of  much  larger  horizontal  extent 
(10  -  >  200  km)  and  greater  intensity  (ion  density  depletions  up  to 
three  orders  of  magnitude)  than  indicated  by  our  small  scale  simulations. 
Therefore,  we  have  extended  our  previous  calculations  and  have  performed 

a  series  of  numerical  simulations  using  much  larger  horizontal  length  scales 

\ 

(~  75km)  for  our  initial  perturbations.  These  seed  long  wavelength 
perturbations,  for  example,  could  be  due  to  neutral  atmosphere  gravity 
wave  effects  [Rottger ,  1976;  Klostermeyer ,  1978;  Booker ,  1979].  At 
the  same  time  we  have  made  very  substantial  improvements  in  the  nu¬ 
merical  techniques  used  to  perform  the  simulations,  including  the 
utilization  of  the  recently  developed  fully  multidimensional  flux- 
corrected  transport  (FCT)  techniques  of  Zalesak  [1979]  .  The  results 
of  our  simulations  indicate  the  following: 

1)  large  horizontal  scale  length  initial  perturbations  evolve 
nonlinearly  into  large  horizontal  scale  length  equatorial 
spread  F  "bubbles;" 

2)  these  bubbles  evolve  on  approximately  the  same  time  scale 
as  do  their  smaller  horizontal  scale  length  counterparts; 


Note:  Manuscript  submitted  December  4,  1979. 
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3)  the  plasma  comprising  these  bubbles  has  its  origin  at  much 
lower  altitudes  than  that  of  the  smaller  horizontal  scale 
length  bubbles,  resulting  in  plasma  density  depletions  of 
very  close  to  100%. 

This  last  result  is  due  to  the  fact  that  the  polarization  fields 
produced  by  ionospheric  irregularities,  aligned  vertically,  possess  a 
fringe  field  component  whose  vertical  extent  is  proportional  to  the 
horizontal  extent  of  the  irregularities  producing  the  field.  This  is 
quite  similar  in  origin  to  the  fringe  field  produced  at  the  edge  of 
a  parallel  plate  capacitor.  Since  the  vertical  extent  of  this  fringe 
field  determines  the  minimum  altitude  from  which  the  rising  bubble  can 
draw  plasma,  it  is  not  surprising  that  larger  horizontal  scale  bubbles 
are  more  severely  depleted.  In  Section  II  we  give  a  brief  review  of  the 
relevant  theory,  and  of  the  basic  equations  used  in  our  simulations. 
Section  III  contains  the  numerical  simulations  and  a  physical  inter¬ 
pretation  of  the  results  is  given.  Section  IV  contains  a  summary,  and 
in  the  appendix  we  describe  briefly  the  numerical  techniques  used  in  our 
present  computer  code,  emphasizing  the  differences  between  the  present 
code  and  our  previous  one. 


II.  THEORY 


The  geometry  of  the  physical  problem  we  are  modeling  is  the  same 
as  in  Ossakow  et  al.  [1979].  All  our  simulations  are  carried  out  in 
a  two  dimensional  (x,y)  coordinate  system.  The  constant  magnetic  field 
_B  is  aligned  along  the  z  axis  (pointing  north).  Gravity  is  directed 
along  the  negative  y  axis.  Since  our  equations  are  two  dimensional, 
n(y),  VR(y),  and  V^(y)  are  taken  to  be  representative  values  of  the 
ambient  electron  density,  recombination  coefficient,  and  ion-neutral 
collision  frequency  (the  result  of  integrating  these  quantities  along 
magnetic  field  lines  and  dividing  by  a  normalizing  scale  length).  Mag¬ 
netic  field  lines  are  assumed  to  terminate  at  both  ends  in  an  electrically 
insulated  medium  (currents  must  close  in  the  two  dimensional  plane,  not 
in  some  distant  E  region). 

Following  Ossakow  et  al.  [1979],  we  describe  the  system  with  the 
two-fluid  plasma  continuity  and  momentum  equations: 
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where  the  subscript  a  denotes  the  species  (i  for  ions,  e  for  electrons), 
n  is  the  species  number  density,  v  is  velocity,  VD  is  the  recombination 
coefficient,  E  is  the  electric  field,  £  is  the  gravitational  acceleration, 
q  is  the  species  charge,  v  is  the  species  collision  frequency  with  the 
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neutral  atmosphere,  U  is  the  neutral  wind  velocity,  c  is  the  speed 

— n 

of  light,  and  m  is  the  species  mass. 

Note  that,  in  contrast  to  Ossakow  et  al  [1979],  we  have  dropped 
the  term  +  VRn^o  from  (1).  This  is  the  equivalent  of  dropping  the 
assumption  of  the  existence  of  an  ionization  source  given  by  that 
term.  This  ionization  source  was  such  that  the  ambient  ionization 
profile  n  (y)  was  an  equilibrium  profile  (dn_^/dt  =  0).  Our  present 
model  therefore  has  instead 


dn 

no 

-  ■  =  -  v_n 

dt  R  ao 
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Hence,  when  normalized  results  n^(x,y) /naQ (y)  are  later  presented, 
it  should  be  understood  that  both  the  numerator  and  denominator  are 
time  dependent. 

Figure  1  shows  the  recombination  rate  VR  and  ion-neutral  collision 

frequency  used  in  our  simulations.  It  shall  be  seen  presently 

that  Vgn  need  not  be  specified  as  long  as  it  is  much  smaller  than 

the  electron  gyro  frequency  ftg.  For  details  on  the  form  of  and 

V  as  depicted  in  Figure  1,  see  Ossakow  et  al.[l979].  If  we  neglect 
R 

the  inertial  terms  (the  left-hand  side)  of  (2)  by  assuming  the  inertial 

time  scales  are  much  larger  than  either  the  gyro  periods  or  the  mean 

time  between  collisions,  then  the  equation  can  be  inverted  to  give  an 

algebraic  expressions  for  v^.  In  two-dimensional  (x,y)  geometry  with 

15  along  the  z  axis,  the  solution  is  for  our  problem,  with  mg  «  , 

V.  Al,  «  1,  V  /ft  =0  (where  ft  =  eB/m  c),  and  U  =0. 
in  i  en  e  a  O  — n 
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The  effect  of  <f>o  is  merely  to  superimpose  a  bulk  westward  plasma 
velocity  on  the  electron  velocity  field  determined  from  $  , 

without  affecting  the  morphology  of  the  developing  structures. 

Hence,  we  ignore  this  motion. 

Our  assumption  of  quasi-neutrality  has  made  one  of  our  two 
continuity  equations  (1)  redundant.  We  therefore  choose  the  electron 
equation  for  its  simplicity: 


(ID 


III.  NUMERICAL  SIMULATION  RESULTS  AND  DISCUSSION 

Equations  (10)  and  (11),  together  with  appropriate  boundary 
conditions,  constitute  the  nonlinear  system  of  equations  we  shall  solve 
numerically.  Note  that  in  contrast  to  Ossakow  et  al.  [1979],  we 
have  chosen  not  to  put  the  equations  into  a  normalized  form  by  dividing 
through  by  nQ(y).  The  numerical  techniques  used  to  solve  these  equa¬ 
tions  are  discussed  in  the  appendix. 

The  numerical  calculations  to  be  presented  were  performed  on 
a  two-dimensional  cartesian  (x,y)  mesh  using  42  points  in  the  x  (east- 
west)  direction,  and  142  points  in  the  y  (vertical) direction.  The 
(uniform)  grid  spacing  was  2km  in  the  y  direction  for  all  calcu¬ 
lations.  The  grid  spacing  in  the  x  direction  was  200m  in  the  "small" 
horizontal  scale  length  cases  and  5km  in  the  "large"  cases.  The 
bottom  of  the  grid  corresponds  to  252km  altitude  and  the  top  of  the 
grid  to  534km  altitude  in  all  simulations.  Periodic  boundary  con- 


ditions  were  imposed  on  both  n  and  <j>^  in  the  x-direction.  In  the  y 
direction  transmittive  boundary  conditions  were  imposed  on  n(dn/dy  =  0) 
and  Neumann  (3<^/dy  =  0)  boundary  conditions  were  imposed  on  <J>^. 

Three  kinds  of  plots  will  be  presented:  (1)  contours  of  con¬ 
stant  n(x,y,t);  (2)  contours  of  constant  n(x,y,t)/n  (y,t);  and  (3)  con- 

o 

tours  of  constant  electrostatic  potential  $  .  Superimposed  on  each 

contour  plot  is  a  dashed  line  depicting  n  (y.t)  for  reference  purposes. 

o 

Our  no(y,o)  profile  is  such  that  the  F2  peak  is  located  at  434km  alti¬ 
tude,  and  the  minimum  electron  density  scale  length  L=n  (dn/dy)  ^ 

■■  o  o 

is  10km.  The  simulations  are  all  identical  except  for  the  east-west 
grid  spacing  Ax  and  the  form  of  the  initial  perturbation.  Two  kinds 
of  initial  perturbations  were  used: 


Perturbation  A: 
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Perturbation  B: 
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Perturbation  A  is  exactly  the  form  used  in  Ossakow  et  al.  [1979], 
and  perturbation  B  is  a  pure  sine  wave  of  wavelength  AO  Ax  (our  system 
length  in  the  x  direction).  Both  represent  maximum  initial  pertur¬ 
bation  amplitudes  of  approximately  5%.  Four  simulations  have  been  run: 
(i)  lS-Perturbation  A  with  Ax  =  200m;  (ii)  lL-Perturbation  A  with 
Ax  =  5km;  (iii)  2S-Perturbation  B  with  Ax  =  200  m;  and  (iv)  2L-Per- 
turbation  B  with  Ax  =  5km.  Calculation  IS  above  is  identical  to 
ESF  III  of  Ossakow  et  al.  [1979].  The  "large"  versus  "small" 
comparison  obviously  involves  comparing  calculation  IS  to  calculation 
1L,  and  calculation  2S  to  calculation  2L.  One  notes  that  for  the 
minimum  L  10km  in  our  simulation,  kL  >  1  for  the  IS  and  2S  cases 
and  kL  <  1  for  the  1L  and  2L  cases. 

Figure  2  shows  isodensity  contours  of  calculation  IS  at  times 
300,  700,  1000,  and  1200  seconds  after  initialization.  Figure  3  shows 
the  same  contours  at  the  same  times  but  for  calculation  1L.  The 
presence  of  much  lower  density  fluid  in  the  bubble  in  calculation  1L 
is  obvious.  Also  obvious  is  a  basic  difference  in  the  bubble  mor¬ 
phology  at  late  times.  At  1200  seconds,  iS  has  pinched  off  into  two 
bubbles,  with  the  more  intense  one  below  the  initial  central  bubble. 

In  addition,  another  bubble  has  formed  in  the  sides  of  the  calcu¬ 
lation  .  These  structures  are  more  obvious  in  the  plot  of  n(x,y)/nQ(y) 
at  1200  seconds  for  IS  shown  in  Fig.  Aa.  The  maximum  depletion  levels 
are  70%  in  the  top  central  bubble  ,97%  in  the  lower  central  bubble, 
and  95%  in  the  side  bubble.  Note  that  here  and  in  all  subsequent 
plots  of  (n,x,z)/nQ(z)  the  contour  plotting  is  such  that  the  first 
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(outer)  depletion  contour  n/n^  is  0.5  and  each  succeeding  inner 

contour  is  0.5  times  the  previous  one.  For  example,  the  lower  bubble 

in  Fig.  4a  has  five  contours.  The  outermost  would  have  n/n  =  0.50 

o 

(50%  depletion),  the  next  inner  one  n/nQ  =  0.25  (75%  depletion),  the 
second  inner  one  n/nQ  =  0.125  (87.5%  depletion)  and  the  innermost 
n/n^  =  0.03  (97%  depletion).  For  the  enhancement  contours  (dashed 
lines)  the  first  outer  contour  is  2.0  and  the  succeeding  inner  ones 
are  2.0  times  the  previous  ones.  In  obtaining  percentage  enhance¬ 
ments  and  depletions  one  then  subtracts  1.0. 

Calculation  1L,  on  the  other  hand,  shows  a  single  plume  of 
depleted  ionization  at  1200  seconds,  with  no  secondary  central 
bubble  and  no  side  bubble.  There  also  is  no  indication  of  a  widen¬ 
ing  of  the  top  of  the  bubble,  as  there  is  in  IS.  In  Fig.  4b  we 
show  a  plot  of  n(x,y)/nQ(y)  for  1L  at  1200  sec.  The  level  of 
depletion  is  greater  than  99.9%  for  the  entire  10km  by  70km  oval 
"hole"  located  inside  the  tenth  solid  contour  of  Fig.  4b  and 
represents  at  least  a  three  order  of  magnitude  decrease  (biteout) 
in  plasma  density. 

We  now  repeat  the  above  comparison,  but  this  time  for  per¬ 
turbation  B  (calculations  2S  and  2L).  Figure  5  shows  isodensity 
contours  of  calculation  2S  at  times  300,  700,  1000,  and  1091 
seconds  after  initialization;  while  Figure  6  shows  similar  plots 
of  calculation  2L  at  times  700,  1000,  1200  and  1364  sec.  Comparison 
again  shows  the  presence  of  much  lower  density  plasma  in  the  bubble 
in  the  2L  calculation.  Morphological  differences  are  also  present, 
the  most  notable  being  the  widening  of  the  top  of  the  bubble  in  2S 


which  is  not  present  in  2L.  Figure  7  shows  a  comparison  of  the  n/n 

o 

profiles  at  late  time.  Again  maximum  depletions  in  the  2S  case  are 
about  97%,  while  a  large  portion  of  the  2L  plume  is  99.9%  depleted  or 
greater. 

We  can  also  compare  the  effect  of  the  form  of  the  perturbation 
by  comparing  IS  to  2S  and  1L  to  2L.  The  latter  comparison  shows  striking 
similarity,  whereas  the  former  shows  some  differences,  the  most  notable 
being  the  lack  of  central  bubble  "pinching"  and  the  lack  of  lateral 
bubbles  in  case  2S.  We  conclude  that  the  morphology  of  the  late-time 
bubbles  is  dependent,  at  least  somewhat,  on  the  form  of  the  initial  per¬ 
turbation. 

Bubble  rise  velocities  are  of  some  interest,  and  we  give  below  the 
average  bubble  rise  velocity  for  each  case,  computed  from  the  last  two 
frames  of  Figs.  2,  3,  5,  and  6: 


IS 

210  m/sec 

1L 

230  m/sec 

2S 

420  m/sec 

2L 

280  m/sec 

The  rise  velocity  of  an  individual  bubble  is  dependent  upon  the  rela¬ 
tive  depletion  level  of  the  bubble,  its  geometry,  and  upon  interactions 
with  other  plasma  structure  nearby.  These  first  two  effects  are  treated 
in  Ossakow  and  Chaturvedi  [1978],  and  the  present  results  above  are  con¬ 
sistent  with  the  results  therein.  For  instance,  the  relatively  high 
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rise  velocity  associated  with  2S  is  seen  to  be  due  to  the  fact  that 


the  bubble  is  more  severely  depleted  than  that  in  IS.  Further,  the 
roughly  equal  rise  velocities  of  the  IS,  1L,  and  2L  bubbles,  in  spite 
of  the  fact  that  the  1L  and  2L  bubbles  are  much  more  severely  depleted 
than  that  in  IS,  is  explained  by  noting  that  IS  actually  approximates 
the  geometry  of  a  "sheet"  of  depleted  plasma,  whereas  the  1L  and  2L 
bubbles  more  closely  resemble  a  cylindrical  geometry  (see  Ossakow 
and  Chaturvedi  [1978]). 

In  an  attempt  to  understand  the  reasons  for  the  differences 
in  the  nonlinear  evolution  of  small  and  large  horizontal  scale  per¬ 
turbations,  we  look  at  the  potential  equation: 


V2 


V(Vinn) 

Vin11 


V  <P  = 


= 

1  C  V 


in 


1  5n 
n  dx 


At  early  times  we  expect  V  <j>  to  be  small  with  respect  to  Bg/cv^, 
so  we  ignore  the  second  term  on  the  left  hand  side,  giving  a  Poisson 
equation  for  (f>^: 

2  A  _  Bg  1  dn 

V  ^1  cV.  n  dx 

in 

We  can  now  interpret  the  right  hand  side  as  simply  the  local  charge 

density  p  (such  that  V.jE  =  p) .  Since  we  have  initialized  all  of  these 
1 

calculations  with  —  —  independent  of  y  (see  (12)),  what  we  are  deal¬ 
ing  with  is  a  distribution  of  charge  density  that  has  the  form  of 
diffuse  "plates"  aligned  in  the  vertical  direction.  Noting  that  the 
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term  v.  decreases  with  altitude,  we  find  that  these  diffuse  "plates" 
in 

have  an  equally  diffuse  "edge"  in  the  y  direction.  Taking  the  analogy 
to  its  conclusion,  we  model  our  initial  conditions,  or  any  vertically 
aligned  structure  for  that  matter,  as  an  array  of  plates  of  charge 
(non-conducting  capacitor  plates)  with  an  edge  somewhere  in  or  above 
our  grid.  Tn  Figure  8  we  show  schematically  the  electric  potential 
field  surrounding  the  edge  of  parallel  plates  of  charge  for  two 
different  separation  distances.  Since  there  is  only  one  scale  length 
in  the  configuration  (the  plate  separation  distance  d),  then  all  other 
scale  lengths  must  be  proportional  to  this  distance.  In  particular, 
the  characteristic  distance  parallel  to  the  plates  over  which  the 
electric  field  outside  the  plates  (the  fringe  field)  falls  off  must  be 
proportional  to  d.  Since  in  our  problem  the  contours  of  electrostatic 
potential  are  in  fact  streamlines  (see  (11)),  this  distance  will 
determine  the  maximum  depth  in  the  fluid  from  which  the  electrostatic 
field  will  draw  fluid  into  the  bubble.  Since  the  plasma  density  is 
lower  at  greater  depths,  this  distance  will  determine  the  minimum 
plasma  density  inside  the  bubble.  To  test  these  ideas,  we  examine 
the  actual  early  time  electrostatic  potential  fields  from  the  cal¬ 
culations  we  have  nresented.  Fig.  9  shows  contours  of  n  and 
for  the  IS  initial  conditions,  and  the  same  plots  for  1L.  A  similar 
comparison  for  cases  2S  and  2L  is  shown  in  Fig.  10.  All  contour 
plots  of  <j>^  are  scaled  in  such  a  way  as  to  evenly  space  exactly  12 
contour  lines  between  the  maximum  and  minimum  value  of  ,  to 
normalize  the  plots  so  they  can  be  compared  without  bias.  The 
increased  vertical  extent  of  the  contours  of  (streamlines)  for 
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cases  1L  and  2L  are  evident. 


Of  course,  the  initial  profile  generating  these  aligned  plates 
of  charge  lasts  only  a  short  time.  In  the  linear  phase  of  growth, 
the  perturbation  grows  in  the  region  of  linear  instability  (the  F  region 
bottomside),  and  damps  elsewhere.  Our  "plates"  therefore  very  quickly 
become  horizontally  spaced  regions  of  charge  with  a  limited  vertical 
extent,  confined  to  the  region  of  steepest  vertical  gradient  on  the  F 
region  bottomside.  Nonetheless,  the  scaling  arguments  advanced  above 
still  hold:  the  vertical  extent  of  the  polarization  electric  field 
scales  as  the  horizontal  extent  of  the  structure  causing  the  field. 

This  is  easily  seen  in  Figure  11  and  12  where  comparison  is  made  of  the 
<j>^  contours  for  IS  vs  1L  and  2S  vs  2L  respectively,  at  a  time  of  700 
sec.  Cases  IS  and  2S  are  seen  to  be  mixing  fluid  over  a  fairly  narrow 
altitude  range,  while  1L  and  2L  have  each  formed  a  large  convective 
cell  more  than  150km  in  vertical  extent,  drawing  plasma  fluid  into  the 
bubble  from  deep  in  the  ionosphere.  It  is  not  surprising,  then,  that 
the  larger  horizontal  scale  bubbles  are  more  severely  depleted  at 
late  times. 

IV.  SUMMARY 

On  the  basis  of  our  numerical  simulations,  and  of  a  qualitative 
scale  analysis  of  the  driving  electrostatic  potential  equation,  we 
conclude  that  the  severe  "biteouts"  of  three  orders  of  magnitude  and 
bubble  rise  velocities  of  150  m/sec  reported  by  McClure  et  al  [1977] 
are  completely  consistent  with  large  ( ~75-200km)  horizontal  bubble 


size  scales.  In  our  simulations,  the  severe  biteouts  associated 
with  large  horizontal  scale  lengths  are  due  to  the  fact  that  the 
plasma  comprising  these  bubbles  has  its  origin  at  much  lower  alti¬ 
tudes  then  in  the  small  horizontal  length  scale  cases.  Again,  these 
results  are  consistent  with  those  of  McClure  et  al.  [1977],  who 
base  their  conclusions  on  ion  mass  spectrometer  measurements  of  the 
H+  -  0+  -  N+  balance  inside  the  bubbles,  which  they  find  to  be 
'^characteristic  of  undisturbed  plasma  found  at  lower  altitudes." 

The  variation  in  the  vertical  velocities  of  various  bubbles  noted 
by  McClure  et  al.  Ll977]  is  probably  due  to  interactions  between 
bubbles.  Note,  for  example,  that  in  Fig.  2,  the  secondary  bubble 
is  rising  at  a  much  slower  rate  than  is  the  central  bubble.  Bubble 
interaction  will  be  the  subject  of  forthcoming  theoretical  and 
numerical  studies. 
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APPENDIX:  Numerical  Solution  of  the  Equations 

Of  the  two  partial  differential  equations  we  must  solve, 

(10)  is  elliptic  and  linear  and  (11)  is  hyperbolic  and  nonlinear. 

Both  equations  represent  numerical  challenges,  and  we  could  easily 
devote  the  bulk  of  this  paper  to  the  numerical  techniques  used  for 
their  solution.  However,  as  we  stated  in  the  introduction,  we  shall 
confine  ourselves  to  the  improvements  made  in  these  techniques  since 
the  calculations  of  Ossakow  et  al.  [1979].  We  begin  with  (11). 

Equation  (11)  is  solved  in  finite  difference  form 
using  a  fully  multidimensional  second  order  in  time,  fourth  order  in 
space,  leapfrog-trapezoidal  flux-corrected  transport  (FCT)  scheme. 

Both  the  higher  order  leapfrog-trapezoidal  scheme  itself,  as  well 
as  the  fully  multidimensional  algorithm  utilized  in  the  critical 
flux-limiting  stage  of  FCT,  are  recent  developments  and  are  described 
by  Zalesak  [1979].  These  developments  represent  significant  extensions 
of  the  theory  of  FCT,  a  numerical  technique  originated  by  Boris  and 
Book  [1973]  to  handle  equations  of  the  form  (11)  where  steep  gradients 
are  expected  for  form.  By  contrast,  the  calculations  in  Ossakow  et  al. 
[1979]  used  a  first  order  in  time,  second  order  in  space  FCT  algorithm 
which  was  only  one-dimensional  (since  fully  multidimensional  FCT 
algorithms  did  not  exist  at  the  time),  and  hence  used  time-splitting 
(sequential  x  and  y  operators)  to  solve  the  two-dimensional  equation 

(11) .  It  is  known  that  time-splitting  can  introduce  numerical 


problems  into  an  incompressible  flow  calculation  (see  Zalesak,  1979), 
although  our  previous  equatorial  spread  F  calculations  did  not 
exhibit  any  of  the  symptoms  of  these  difficulties. 

Equation  (10)  is  the  elliptic  equation  whose  solution 
gives  us  the  electrostatic  potential  $  .  The  right  hand  side  is 
known  and  the  left  hand  side  represents  a  Hermitian  operator 
operating  on  <f>^,  giving  only  real  eigenvalues  and  apparently  easing 
the  difficulty  of  solution.  However,  the  extremely  large  range  for 
the  values  of  the  quantity  n  makes  for  an  equally  large  span  of 
operator  eigenvalues,  and  solution  of  the  equation  in  this  form 
using  iterative  techniques  has  not  been  successful.  We  have  found 
one  (and  only  one)  method  of  direct  solution,  the  stabilized  error 
vector  propagation  (SEVP)  method  of  Madala  [1978],  but  the  execution 
speeds  for  SEVP  are  not  as  favorable  as  for  the  method  we  now  describe. 

We  start  by  expanding  the  operator  and  dividing  through 
by  v.  n,  as  was  done  by  Ossakow  et  al.  [1979],  giving 


v(v.  n) 
in 


V*  =  -  I 

a  1  c  V.  n  dx 
in 


(A-l) 


The  equation  is  now  in  a  form  suitable  for  solution  by  the  Chebychev- 
iterative  relaxation  technique  of  McDonald  [1977].  However,  great 

care  must  be  given  to  the  differencing  of  the  term  ~  -  V(v^nn)  and 

1  ^  1^ 

—  -r—  ,and  this  is  the  point  we  wish  to  address.  We  work  with  the 
n  dx 

1  dn 

term  —  in  one  spatial  dimension,  since  this  example  will  make 

the  point.  A  straightforward  second-order  difference  form  for  this 
term  is 
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1 


a  =  j_  ni+rni-i 
'n  3x  *  2Ax  n^ 


(A-2) 


where  the  subscript  i  refers  to  grid  point  location  in  the  x  direction 
and  Ax  is  the  (uniform)  spacing  between  grid  points.  This  is  the 
form  used  in  Ossakow  et  al.  [1979].  We  shall  show  that  this  difference 
form  produces  solutions  <j>^  with  potentially  undesirable  properties, 

and  causes  undue  numerical  hardship  in  finding  those  solutions. 

1  dn  B 

Let  us  rewrite  —  as  -r—  (In  n).  A  physical  inter- 

n  ox  ox 

pretation  of  the  term  is  now  much  easier:  the  term  tells  us  how 

rapidly  the  logarithm  of  n  is  varying  with  respect  to  x.  Suppose, 

for  argument's  sake,  that  the  smallest  and  largest  values  of  n  in 

the  problem  are  10^  and  10^  respectively.  On  a  grid  of  size  Ax',  the 

largest  value  representable  for  that  term  would  occur  when  a  fluid 

element  of  density  101  and  one  of  density  105  occupy  adjacent  grid 

£ 

points.  The  value  of  —  (In  n)  evaluated  midway  between  these  two 

ox 

grid  points  would  be  ~  (  In  10^  -  In  101)  =  9.2  Ax  Evaluation 
of  (A-2)  for  n^  ^,n^,  an<*  ni+l  ^av^n8  values  of  10  ,  10^,  and  10^ 
respectively  gives  a  value  for  (^  -^)  of  5  x  10^  Ax  \  far  in 
excess  of  the  maximum  value  for  this  term  given  by  the  above  argu¬ 
ment.  Logarithmic  interpretation  of  this  term  would  state  that  n 

4 

varied  by  more  than  10  orders  of  magnitude  over  a  single  grid 
spacing,  a  ridiculous  statement  in  light  of  the  fact  that  there 
are  only  four  orders  of  magnitude  of  n  in  the  problem. 
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of  the  term 


The  potentially  large  values  of  these  terms,  in  particular 


-  V  (v,  n)  in  (A-l),  not  only  cause  extremely  slow 

V,  n  j.  in  } 

in 

convergence  of  the  iterative  solution,  but  can  also  put  spurious 


oscillations  into  the  exact  finite  difference  solution  ^  due  to 
cell  Reynolds  number  effects  [Roache,  1976].  As  shown  by  Roache 


[1976]  these  oscillations  can  occur  any  time  the  value  of  the 


term  -  V  (V.  n)  in  (A-l)  exceeds  a  critical  value  of  2  Ax 

V.  n  jl  in 
in 

It  is  clear,  then,  that  (A-2)  represents  an  undesirable  difference 


form  for  these  logarithmic  terms.  Better  approximations  are 


(I  ^  .  _L_ 

'n  dz  ^  2Ax 


(In  n  -  In  n  ) 


(A- 3) 


,1  dnN  1  "i+l  "  "i-l 

and  (—  =  -t—  - - - 

n  dx  t  Ax  n1+1  +  n  ^ 


(A-4) 


Equation  (A-3)  is  probably  the  most  accurate,  but  evaluation  of 
the  logarithms  at  every  time  step  is  computationally  expensive,  and 
for  problems  like  the  one  at  hand  where  n  varies  by  many  orders  of 
magnitude  across  the  grid,  there  is  still  no  guarantee  that  the 
critical  cell  Reynolds  number  will  not  be  exceeded.  For  these 
reasons  we  use  (A-4)  for  the  present  calculations. 
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ION-NEUTRAL  COLLISION  FREQUENCY,  ^(secf1) 

e  ion-neutral  collision  frequency  Vin  and  recombination  rate 
function  of  altitude)  used  in  the  numerical  simulations 


Fig.  2  -  Sequence  of  four  plots  showing  iso-electron  density 
contours  of  calculation  IS  at  300,  700,  1000,  and  1200  sec. 
Superimposed  on  each  plot  is  a  dashed  line  depicting  n  (z,t). 
Electron  densities  are  given  in  cm  .  The  observer  is  looking 
southward . 
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Same  as  Fig.  2  but  for  calculation  2S  at  300,  700 
1000,  and  1091  sec. 
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Fig.  6  -  Same  as  Fig.  5  but  for  calculation  2L  at  700,  1000, 

1200,  and  1364  sec. 
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Schematic  representation  of  contours  of  constant  electrostatic  potential  produced  at 
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trie  field  found  between  the  plates.  This  facilitates  the  comparison  of  the  rate  at  which  the 
electric  field  falls  off  away  from  the  edge  of  the  plates.  The  much  more  rapid  falloff  in  a) 
is  evident. 
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Fig.  10  -  Same  as  Fig.  9  but  for  cases  2S  and  2L 
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OICY  ATTN  COOE  R720  J.  WORTHINGTON 


DIRECTOR 

JOINT  ST RAT  TGT  PLAAMlNG  STAFF 
OFFUTT  AFB 
OMAHA,  NB  68113 

OICY  ATTN  JLTW-2 
OICY  ATTN  JP$T  G.  GOETZ 

CHIEF 

LIVERMORE  DIVISION  FLO  COmAAO  DMA 
DEPARTMENT  OF  DEFENSE 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94SS0 
OICY  ATTN  FCPRl 


01 RECTOR 

NATIONAL  SECURITY  AGENCY 

DEPARTMENT  OF  DEFENSE 

FT.  GEORGE  G.  MEADE,  MD  207SS 

OICY  ATTN  JOHN  SKILLMAN  R52 
OICY  ATTN  FRANC  LEONARD 
OICY  ATTN  Wl4  PAT  CLARK 
OICY  ATTN  OLIVER  H.  BARTLETT  W32 
OICY  ATTN  R5 

COMMANDANT 

NATO  SCHOOL  (SHAPE) 

APO  NEW  YORK  09172 

OICY  ATTN  U.S.  DOCUMENTS  OFFICER 

UNDER  SECY  OF  DEF  FOR  RSCH  8  ENGRG 
DEPARTMENT  OF  DEFENSE 
WASHINGTON,  O.C.  20301 

OICY  ATTN  STRATEGIC  6  SPACE  SYSTEMS  (OS) 


DIRECTOR 

DEFENSE  COMMUNICATIONS  AGENCY 
WASHINGTON,  O.C.  20305 

(AOR  CNWDI:  ATTN  COOE  240  FOR) 


OICY  ATTN  COOE  1018 


DEFENSE  DOCUMENTATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA.  22314 

(12  COPIES  IF  OPEN  PUBLICATION,  OTHERWISE  2  COPIES) 
12CY  ATTN  TC 


DIRECTOR 

DEFENSE  INTELLIGENCE  AGENCY 
WASHINGTON,  D.C.  20301 
OICY  ATTN  DT-1B 
OICY  ATTN  D0-4C  E.  O' FARRELL 
OICY  ATTN  DIA AP  A.  WISE 
OICY  ATTN  DIAST-5 
OICY  ATTN  DT-1BZ  R.  MORTON 
OICY  ATTN  HQ-TR  J.  STEWART 
OICY  ATTN  W.  WITTIG  OC-7D 

DIRECTOR 

DEFENSE  NUCLEAR  AGENCY 
WASHINGTON,  D.C.  20305 
OICY  ATTN  STVL 
04CY  ATTN  T|TL 
OICY  ATTN  DOST 
03CY  ATTN  RAAE 

COMMAFOER 
FIELD  COftlANO 
DEFENSE  NUCLEAR  AGENCY 
KIRTLAND  AF3,  NM  87115 
OICY  ATTN  FCPR 

DIRECTOR 

INTERSERVICE  NUCLEAR  WEAPONS  SCHOOL 
KIRTLAND  AFB,  *1  87115 

OICY  ATTN  DOCUMENT  CONTROL 


WWMCCS  system  engineering  org 
WASHINGTON,  D.C.  20305 

OICY  ATTN  R.  CRAWFORD 

COMMAFOER /DIRECTOR 
ATMOSPHERIC  SCIENCES  LABORATORY 
U.S.  ARMY  ELECTRONICS  COMMAND 
WHITE  SA*)S  MISSILE  RANGE,  NM  88002 
OICY  ATTN  DELAS-EO  F.  NILES 

DIRECTOR 

BFf)  ADVANCED  TECH  CTR 
HUNTSVILLE  OFFICE 
P.  0.  BOX  1500 
HUNTSVILLE,  AL  35807 

OICY  ATTN  ATC-T  MELVIN  T.  CAPPS 
OICY  ATTN  ATC-0  W.  DAVIES 
OICY  ATTN  ATC-R  DON  RUSS 

PROGRAM  MANAGER 
BMD  PROGRAM  OFFICE 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22333 

OICY  ATTN  DACS-3MT  J.  SHEA 

CHIEF  C-E  SERVICES  DIVISION 
U.S.  ARMY  COMMUNICATIONS  CFO 
PENTAGON  RM  1B269 
WASHINGTON,  D.C.  20310 

OICY  ATTN  C-E-SERVICES  DIVISION 


COMMANDER 

FRADCOM  TECHNICAL  SUPPORT  ACTIVITY 
DEPARTMENT  OF  THE  ARMY 
FORT  MONMOUTH,  N.J.  07703 

OICY  ATTN  DRSEL-NL-RD  H.  BENNET 
OICY  ATTN  DRScL-PL-ENV  H.  0OMCE 
OICY  ATTN  J.  E.  QUIGLEY 


JOINT  CHIEFS  OF  STAFF 
WASHINGTON,  D.C.  20301 

OICY  ATTN  J-3  WWMCCS  EVALUATION  OFFICE 
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COMMANDER 

HARRY  DIAMOND  LABORATORIES 
DEPARTFCNT  OF  THE  ARMY 
2800  POWDER  MILL  ROAD 

AOELPMI,  MO  20785  ,  _ 

(CNWDI-1W€R  envelope:  ATTN:  OELK>-RBh) 
01CY  ATTN  OELtO-TI  M.  WEINER 

OICY  ATTN  OELTC-RB  R.  WILLIAMS 

01CY  ATTN  OELTO-NP  F.  WlXCNITZ 

01CY  ATTN  OELTO-NP  C.  MOAZED 


COMMANDER 

U.S.  ARMY  COTTt-ELEC  ENGRG  INSTAL  ASY 
FT,  rlUACHUCA,  A2  8  5615 

01CY  ATTN  CCC-EMEO  GEORGE  LAX£ 


COMMANDER 

J.S.  ARMY  FOREIGN  SCIENCE  (  TECH  CTR 
220  7TM  STREET,  NE 
Charlottesville,  va  22901 

0  ICY  ATTN  DRXST-SO 
01CY  ATTN  R.  JONCS 


COMMANDER 

U.S.  ARMY  MATERIEL  DEV  8  READINESS  CNC 
5001  EISENHOWER  AVENUE 
ALEXANDRIA,  VA  22555 

01CY  ATTN  DRCLDC  J.  A.  BENDER 


CONMANOER 

L.s.  ARMY  NUCLEAR  AND  CHEMICAL  AGENCY 
7500  BACKLICK  ROAD 
BLOG  2075 

SPRINGFIELD,  VA  22150 
01CY  ATTN  LIBRARY 

DIRECTOR 

U.S.  ARMY  BALLISTIC  RESEARCH  LABS 
ABERDEEN  PROVING  GROUND,  MD  21005 

01CY  ATTN  TECH  LIB  EDWARD  BA ICY 

CONMANCER 

U.S.  ARMY  SATCOM  AGENCY 
FT.  MONFDUTH,  NJ  07705 

01CY  ATTN  DOCUMENT  CONTROL 

COMMANDER 

U.S.  ARMY  MISSILE  INTELLIGENCE  AGENCY 
REDSTONE  ARSENAL,  At  55809 
OICY  ATTN  JIM  GAMBLE 


U.S.  army  traooc  systems  analysis  activity 
WHITE  SANDS  MISSILE  RANGE,  NM  88002 
OICY  ATTN  ATAA-SA 
OICY  ATTN  TCC/F .  PAYAN  UR. 

OICY  ATTN  ATAA-TAC  LTC  U.  HESSE 

COMMANDER 

NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
WASHINGTON,  D.C.  20560 

OICY  ATTN  NAVALEX  054  T.  HUGHES 
0 ICT  ATTN  PME  117 
OICY  ATTN  PME  117-T 
OICY  ATTN  CODE  5011 

COMMANDING  OFFICER 
NAVAL  INTELLIGENCE  SUPPORT  CTR 
4501  SUITLANO  ROAD,  BLDG.  5 
WASHINGTON,  D.C.  20390 

OICY  ATTN  MR.  DUBBIN  STIC  12 

OICY  ATTN  NlSC-50 

OICY  ATTN  CODE  5404  U.  GALET 


COMMANDER 

naval  surface  weapons  center 

DAH.GREN  LABORATORY 
DAHLGREN,  VA  22448 

OICY  ATTN  CODE  OF -14  R.  SUTLER 

COMMANDING  OFFICER 
NAVY  SPACE  SYSTEMS  ACTIVITY 
P.O.  BOX  92960 
WORLOWAY  POSTAL  CENTER 
LOS  ANGELES,  CA.  90009 
OICY  ATTN  CODE  52 

OFFICE  OF  NAVAL  RESEARCH 
ARLINGTON,  VA  22217 

01CT  ATTN  CODE  465 
OICY  ATTN  CODE  461 
01CT  ATTN  CODE  402 
OICY  ATTN  CODE  420 
OICY  ATTN  COOE  421 

CONMMDER 

AEROSPACE  DEFENSE  COHMAND/DC 
DEPARTNtNT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

OICY  ATTN  DC  MR.  LONG 

COMMANDER 

AEROSPACE  DEFENSE  COMMANO/XPD 
DEPARTMENT  OF  THE  AIR  FORCE 
ENT  AFB,  CO  80912 

01CT  ATTN  XPOQQ 

OICY  ATTN  XP 

AIR  FORCE  GEOPHYSICS  LABORATORY 
HANSC0M  AFB,  MA  01751 

OICY  ATTN  OPR  HAROLD  GARDNER 

OICY  ATTN  OPR-1  JAMES  C.  ULW1CX 

01CT  ATTN  LKB  KEWVETH  S.  W.  CHAMPION 

OICY  ATTN  OPR  ALVA  T.  STAIR 

OICY  ATTN  PHP  JULES  AARONS 

OICY  ATTN  Pro  JURGEN  BUCHAU 

OICY  ATTN  PNC  JOFN  P.  MULLEN 

AF  WEAPONS  LABORATORY 
XIRTLAND  AFB,  NM  87117 
OICY  ATTN  SUL 

OICY  ATTN  CA  ARTHUR  H.  GUENTHER 
OICY  ATTN  OTC  CAPT  J.  BARRY 
OICY  ATTN  DYC  JOHN  M.  KANM 
OICY  ATTN  DYT  CAPT  MARK  A.  FRY 
OICY  ATTN  DES  MAJ  GARY  GANONG 
OICY  ATTN  DYC  J.  JANNI 

AFTAC 

PATRICK  AFB,  FL  52925 

OICY  ATTN  TF/MAJ  WILEY 
OICY  ATTN  TN 

AIR  FORCE  AVIONICS  LABORATORY 
WRIGHT-PATTERSON  AFB,  OH  45455 
OICY  ATTN  AAD  WADE  HUNT 
OICY  ATTN  AAD  ALLEN  JOHNSON 

DEPUTY  CHIEF  OF  STAFF 
RESEARCH,  DEVELOPMENT,  8  ACQ 
DEPARTMENT  OF  THE  AIR  FORCE 
WASHINGTON,  D.C.  20550 
OICY  ATTN  AFRDQ 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIV1SION/XR 
DEPARTMENT  OP  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01751 

OICY  ATTN  XR  J.  DEAS 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIV1S10N/YSEA 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01751 
OICY  ATTN  YSEA 


COMMANDED 

NAVAL  ocean  systems  center 
SAN  0 I C^O,  CA  92152 

OKt  ATTN  CODE  552  *.  MGlE* 

OICY  ATTN  CODE  0230  C.  8AGGETT 
OICY  ATTN  CODE  6  1  R.  EASTON 

01  RECTOR 

NAVAL  RESEARCH  LABORATORY 
wAShINGtON,  J.C.  20375 

01CY  ATTN  COOE  6700  TIMOTHY  P.  COFFEY 

V.2S  CYS  IP  UNC-ASS,  1  CY  IF  CLASS ) 

OICy  ATTN  COOE  6701  JACK  D.  BRCWN 
01CY  ATTN  CODE  67110  BRANCH  HEAD  Cl 50  CYS 
JP  oNCLASS,  1  CY  IF  CLASS) 

*1CY  ATTN  CODE  7500  hQ  COMM  DIR  BRUCE  wAlD 

01CY  ATTN  COOE  7550  J.  DAVIS 

OICY  ATTN  CODE  7580 

01CY  ATTN  COOE  7551 

OICY  ATTN  COOE  7555 

OICY  ATTN  CODE  6730  E.  MClEAn 

OiCY  ATTN  COOE  7127  C.  JOHNSON 

COMMANDER 

NAVAL  sea  systems  command 
WASHINGTON,  O.C.  20362 

OICY  ATTN  CAPT  R.  PITKIN 

COMMANDER 

naval  space  surveillance  system 

OAHLGREN,  VA  22448 

OICY  ATTN  CAPT  J.  H.  BuRTON 

OFFICER -IN-CHARGE 

naval  surface  weapons  center 
white  OAK,  SILVER  SPRING,  MO  20910 
OICY  ATTN  CODE  FJi 

DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  OFFICE 
DEPARTMENT  of  The  navy 
WASHINGTON,  O.C.  20376 
OICY  ATTN  NSP-2141 
C\CY  ATTN  NSSP-2722  fREO  -WIMBERLY 

NAVAL  SPACE  SYSTEM  ACTIVITY 
P.  O.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALIF.  90009 

G1CV  ATTN  A.  8.  HAZZARO 

HEADQUARTERS 

ELECTRONIC  SYSTEMS  DIVISION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AF3,  MA  01731 

OiCY  ATTN  DOCC  MAJ  J.  C.  CLARK 


COMMA  NUck 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WRIGHT-PATTERSON  AF8,  OH  45433 
OiCY  ATTN  NICD  LIBRARY 
OICY  ATTN  ETDP  B .  BALLARD 

COMMA  I'D  ER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
GRIFF ISS  AF0,  NY  13441 

OICY  ATTN  DOC  LIBRARY /TSLD 
OICY  ATTN  UCSE  V.  COYNE 

FAMSO/SZ 

POST  OFFICE  BOX  92960 
WO«L>AY  POSTAL  CENTER 
LOS  ANGclES,  CA  90009 
C SPACE  DEFENSE  5vST£MS) 

OICY  ATTN  SZj 


STRATEGIC  AM  COMMAND/ *PF S 
OFFuTT  AFB,  NB  68113 

OICY  ATTN  XPFS  MAj  3.  STEPHAN 
OICY  ATTN  ADwATE  MAj  BRv/CE  BAuER 
OICY  ATTN  NRT 

OICY  ATTN  DOK  CHIEF  SCIENTIST 

SAMSO/YA 

P.  O.  BOX  92960 

WORLDWAY  POSTAL  CENTER 

.-OS  ANGELES,  CA  90009 

OICY  ATTN  CAPT  L.  BLACKwELDER 

SAMSO/SK 
p.  0.  BOX  92960 
worldway  postal  center 
LOS  ANGELES,  CA  90009 

OICY  ATTN  SKA  (SPACE  COMM  SYSTEMS)  M.  CLAVIN 
SAMSO/MN 

NORTON  AFB,  CA  92409 
(M1NUTEMAN) 

OICY  ATTN  MNNL  lTC  KENNEDY 
COMMANDER 

ROME  AIR  DEVELOPMENT  CENTER,  AFSC 
HANSCOM  AFB,  M&  01731 

OICY  ATTN  EEP  A.  LORENTZEN 


DEPARTMENT  OF  ENERGY 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  5400 
Al3uQuERQuE,  NM  67115 

OICY  ATTN  DOC  CON  FOR  D.  ShERwOOO 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C .  20545 

OICY  ATTN  DOC  CON  FOR  A.  lABOwITZ 

EG6G,  INC. 

LOS  ALAMOS  DIVISION 

p.  0.  BOX  809 

LOS  alamos,  nm  85544 

OICY  ATTN  DOC  CON  FOR  J.  BPEEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORv 
p.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  T£Ch  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  u-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  u-31  R.  HAGER 
OICY  ATTN  DOC  CON  FQR  t-46  F.  SEWARD 

LOS  alamos  scientific  laboratory 
P.  0.  60X  1663 

los  alamos,  nm  87545 

OICY  ATTN  DOC  CON  FOR  J.  WOLCOTT 
OICY  ATTN  DOC  CON  FOR  R.  F.  TASCHEK 
OICY  ATTN  DOC  CON  FOR  E.  JONES 
OICY  ATTN  DOC  CON  FOR  j.  MALlK 
OICY  ATTN  DOC  CON  FOR  R.  JEFFRIES 
OICY  ATTN  DOC  CON  FOR  J.  ZINN 
OICY  ATTN  DOC  CON  FOR  P.  KEATON 
OiCY  ATTN  DOC  CON  FOR  D.  wESTERVElT 

SAND I A  laboratories 

p.  0.  BOX  5800 
ALBUQUERQUE ,  NM  87115 

OICY  ATTN  DOC  CON  FOR  J.  MARTIN 

OICY  ATTN  DOC  CON  FOR  W.  BROWN 

OICY  ATTN  DOC  CON  FOR  A.  ThCRNBROUGh 

QICY  ATTN  DOC  CON  FOR  T.  wRIGhT 

OICY  ATTN  DOC  CON  FOR  D.  OAhLGREN 

OICY  ATTN  DOC  ON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 
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SAfV.UA.  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  dOX  969 
LIVERMORE,  CA  94550 

01CY  ATT\  JOC  CON  FOR  d.  MuRPhEy 
OICY  ATTn  JOC  CON  FOR  T.  COC* 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
*A$nlNGTuN,  D.C.  20545 

01CY  ATTN  JOC  CON  FOR  D.  GAl£ 


O'mER  GOVERNMENT 


central  INTELLIGENCE  AGENCY 
ATTN  RD/SI#  RM  5G48,  hQ  BLOG 
WASHINGTON,  O.C.  20505 

01CY  ATTN  OSI/PSIO  RM  5P  19 

DEPARTMENT  of  COMMERCE 
national  Bureau  of  standards 
WASHINGTON,  O.C.  20254 

Call  corres:  attn  sec  officer  for) 

01CY  ATTN  R.  MOORE 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  S  INFO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CR0M6IE 
01CY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  $  ATMOSPHERIC  ADMIN 

environmental  research  laboratories 

DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 

OICY  mTTN  R.  GRUBB 

01CY  ATTN  AERONOMY  LAB  G.  REID 


AEROSPACE  CORPORATION 

P.  0.  BOX  92957 

uOS  ANGELES,  CA  9U0Q9 

OICY  ATTN  I.  GARFUISKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  N.  STOCKWELL 
OICY  ATTN  D.  OLSEN 
OICY  ATTN  J.  CARTER 
JlCY  ATTN  F .  MORSE 
OICY  ATTN  SMFA  FOR  PWW 

ANALYTICAL  SYSTEMS  ENGI BERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  ilA  01803 

OICY  ATTN  RADIO  SCIENCES 

BERKLEY  RESEARCH  ASSOCIATES,  INC. 
P.  0.  BOX  983 
BERKELEY,  CA  94701 

OICY  ATTN  J.  WORKMAN 


BOEING  COMPANY,  TnE 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

OICY  ATTN  G.  KEISTER 
OICY  ATTN  o.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KEINMEY 

CALIFORNIA  AT  SAN  DIEGO,  UNI  V  OF 
1PAPS,  B-019 
LA  JOLLA,  CA  92093 

OICY  ATTN  HENRY  G.  BOCXER 

BROWN  ENGINEERING  COMPANY ,  INC. 

Cummings  research  park 
HUNTSVILLE,  AL  35807 

OICY  attn  ROMEO  A.  DtLlBERlS 

Charles  stark  draper  laboratory,  inc. 

555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 

OICY  ATTN  D.  8.  COX 
OICY  ATTN  J.  P.  GILMORE 

CO*>uTER  SCIENCES  CORPORATION 
6565  ARLINGTON  3LVD 
FALLS  CHuRCn,  VA  2 2046 
oicy  attn  h .  blank 
OICY  ATTN  JOHN  SPOOR 
OICY  ATTN  C.  NAIL 

COMSAT  LABORATORIES 
L I  NTH 1 CUM  ROAD 
CLARKSBURG,  MD  20734 

OICY  ATTN  G.  HYDE 

CORNELL  JNIVERSITY 

department  of  electrical  engineering 
ithaca,  NY  14850 

OICY  ATTN  D.  T.  FARLEY  JR 

electrospace  systems,  inc. 
box  1359 

RICHARDSON,  TX  75080 

OICY  ATTN  h.  logston 

OICY  ATTN  SECURITY  (.PAuL  PHILLIPS) 

ESl  INC. 

495  JAVA  DRIVE 
Sunnyvale,  CA  94086 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
OICY  ATTN  C.  W.  PRETTIE 

FORD  AEROSPACE  (  COMMUNICATIONS  CORP 
3939  FABIAN  WAY 

palO  alto,  ca  94303 

OICY  ATTN  J.  T.  MATT 1NGLEY 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  space  center 
GOOOARD  BLVO  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  SORTER  SPACE  SCI  lA0 

GENERAL  ELECTRIC  CO*>ANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  1J201 

OICY  ATTN  F.  RE  I  BERT 

general  electric  cok>any 
te*>o-center  for  aovanceo  stuoies 

816  STATE  STREET  (P.O.  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OICY  ATTN  DASIAC 

oicy  atyn  don  chandler 

OICY  ATTN  TDM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  KNA>P 
OICY  ATTN  WILLIAM  MCNAMARA 
C1CY  ATTN  8.  GAMBlLL 
OICY  ATTN  MACK  ST  ANT  > 
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J  *c. 


' r  ScRtlLcS  CD.  , 

-.vcS 

C.XRT  $:*££’ 
s'rR-C.Sc,  \r  1  5201 

)!-■*  AT'\  u.  «i  w_  V\ 

i,E\ERAi_  REScARCm  CORPORATION 
SANTA  BARBARA  DIVISION 
p.  o.  box  6770 
SANTA  dARSARA,  CA  93111 

01 CY  ATTN  wOnb  ISE  JR 
OICY  ATTN  JOEu  GAR6ARINO 

geophysical  institute 
university  of  Alaska 

FAIRBANKS,  AK  99701 

(ALk.  Class  attn:  security  officer) 

OICY  ATTN  T.  N.  DAVIS  (uNCL  ONLY) 

OICY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

electronics  systems  grp-eastern  oiv 

77  A  STREET 
NEEDHAM,  MA  02194 

OICY  ATTN  MARSHAL  CROSS 

Illinois,  university  of 
department  of  electrical  engi^ering 

URBANA,  :l  61803 

OICY  ATTN  K.  TEH 

Illinois,  university  of 

107  COBLE  HALL 
801  S.  WRIGHT  STREET 
URBANA,  1l  60680 

(All  CORRES  ATTN  SECURITY  SUPERVISOR  FOR) 
OICY  ATTN  K.  YEH 


KAMAN  SCIENCES  COR3 
■*.  0.  BOx  7 -*63 
COLORADO  SPRINGS,  CD  60953 
OICY  ATTN  T.  MEAGHER 

-1  MCA* IT  COR? 

10453  ROSELLE 

SAN  D1ESO,  CA  92121 

OICY  ATTN  19*1 \  JACOBS 


M.I.T.  LINCOLN  LABORATORY 
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